function [dat]=r_mld
% ---- Read BATS parameter

% ---- Create the Structure ...        
dat = struct('dummy',1);

dirMLD='/calcul/pc-charria/OBSERVATIONS/CLIMATO/MLD/';

ncload([dirMLD,'mldR003_clim_1967_2006.nc'],'lon','lat','mld')
dat.londata=lon;
dat.latdata=lat;

dat.londata(find(dat.londata > 180))=dat.londata(find(dat.londata > 180))-360;
[nn,I]=sort(dat.londata);
dat.londata=nn;
dat.data=mld(:,:,I);
dat.data=permute(dat.data,[3 2 1]);
dat.timedata=1:12;

% ---- Clean the Structure ...
dat=rmfield(dat,'dummy');    

